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Abstract. Let / be a classical holomorphic cusp form for SL^^E) of weight k which is a normalized cigen- 
function for the Hecke algebra, and let A(n) be its eigenvalues. In this paper we study "shifted convolution 
sums" \(n)\(n + h) after averaging over many shifts h and obtain asymptotic estimates. The result is 
somewhat surprising: one encounters a transition region depending on the ratio of the square of the length 
of the average over h to the length of the shifted convolution sum. The phenomenon is similar to that 
encountered by Conrey, Farmer and Soundararajan in their 2000 paper pQ on transition mean values of 
the Jacobi symbol, and the connection of both results to Eisenstein series and multiple Dirichlet series is 
discussed. 

Let / be a classical holomorphic Hcckc eigencuspform for SL-itffc) of weight k G 2Z, k > 0. It admits a 
Fourier expansion of the form 

f(z) = ^n^A/(n)eH, 

n>l 

where the A/(n) are the Hecke eigenvalues of /, normalized so that A/(l) = 1, and we use the standard 
notation e(z) := e 27riz . In this paper, we study "shifted convolution sums", i.e. sums of the form 

5^A/(n)A/(n + /i). 

n 

These sums have many applications in analytic number theory. They often arise in subconvexity results for 
GL\ and GL2 L-functions, see Sarnak [14], or Lecture 4 of Michel in [13], and moreover, they are crucial 
to Soundararajan and Holowinsky's resolution of the quantum unique crgodicity conjecture for classical 
holomorphic modular forms, see [9]. In this paper we obtain asymptotic estimates for shifted convolution 
sums after averaging over many shifts h, i.e. we study sums of the form 

(1) £ E A /(»)M» + *0, 

where the notation n x X indicates a sum over n of length X with a for-now-unspecified smoothing. 

We show that when X and Y grow large in such a manner that Y 2 /X 00, the double sum (JT|) has 
an asymptotic formula with well controlled error terms, and we obtain nontrivial asymptotic upper bounds 
when Y 2 /X — > 0. However, the most interesting case of our result is the transition phase between these 
two situations, that is, when X and Y go to infinity and Y 2 jX = c, a fixed constant. In this situation, the 
asymptotic we prove depends delicately on the constant c. One perspective is to interpret the sum (Q]) as 
varying on the open first quadrant of the (X, y 2 )-plane. Asymptotic estimates as X and Y 2 go to infinity 
are a description of the singularity at infinity in this quarter-plane. In this paper we find that the asymptotic 
behavior varies continuously on a blowup of the point at infinity in the quarter-(A", y 2 )-plane. Our results 
for shifted convolution sums are very similar to the interesting results of Conrey, Farmer and Soundararajan 
on transition mean values of the Jacobi symbol Q] . They study the sum 



S(X,Y):= £ £ ( 



%<X n<Y 
1 odd n odd 



/TO 

V n 



and similarly find asymptotic formulae when one of either X or Y grow much faster than the other, and a 
transition region when X/Y is a fixed constant, in which the asymptotic varies continuously on a blowup. 
The sums ^ and S(X,Y) at first look dissimilar. Note however that the problem of determining the 



asymptotic behavior of ([!} is equivalent to determining asymptotic behavior in the sum 

2 



E ( E A/(n + m) 



which puts ([1} and S(X, Y) on equal footing, and makes apparent why the transition region is at Y 2 / X = c. 
Moreover, we will see later that these two examples rest on the same ideas and have many common features. 
First, however, we state precisely a corollary of the Main Theorem (found in section 4) of this paper. The 
Main Theorem is very similar to the corollary, but allows a choice of cut-off functions. 
We define the function cj on the positive real numbers 



c /0) : = ~y a Yl x f( n ) 2w k {n 2 na) 



n>l 

with 



1 f T(s + k - 1)T (s - k) 

W k (x) := — / -i — ±-£ 2 ->-x- s ds 

2 ™ J(i+a) r ( 2 - s) 

for any fixed a > 0. 

Corollary 1. If 1 <Y < X then 

£ S>MM« + *> (^f e r^ - («, (£) - m^H y + o t (xiYiW) 

where Cf(a) defined above is a smooth function on the positive real numbers which 

• as a — > oo decays faster than any polynomial 

• as a — > is 

r(fc)L(l,sym 2 /) 
= 2C(2l + Sfe(a) ' 

where Ek(a) = ( a ^j unconditionally, and Ek{a) = Ok.e (ai~ E ^j for any e > assuming the 

Riemann hypothesis for the classical Riemann zeta function, 
and where 9 = or = ^ depending on whether one assumes the generalized Ramanujan conjecture for Maass 
forms of SL2 (Z) or not. 

As a benchmark, the best point-wise estimates for shifted convolution sums on SL2 (Z) give 

(2) Yl *f(n)\f(n + h)< e X*+*, 

see Sarnak [14] . In Corollary 1, observe that if Y 2 is large compared to X we have that 
\- , ( u , , 1 \ f njn + h) ^ ^ r(fc)L(l,sym 2 /) y 

h<Yn>l v 7 SK ' 

and if Y 2 is small compared to X then 

E E A/(n)A /( n + *) f^T) ^ = Cfc (jtfy) + , 

d<y»>i ^ 7 

or an even better bound if we assume the Riemann Hypothesis. In the transition region when Y 2 is a 
constant multiple of X, the asymptotic growth is controlled by the function c/(a). 
Now we describe the work of Conrey, Farmer and Soundararajan. Their result is 



Theorem (Conrey, Farmer and Soundararajan). Uniformly for all large X and Y , we have 



S(X,Y) = -^C ( ] +0 [(XYts +YX^)logXY) , 



f2nk 2 \ . f2irk 2 \ , 
cos + sin I ay. 



where for a > we define 

1 00 1 pa / 

j4n alternate expression for C(a) is 

C(a) = a + a* I g 1 ^ ^sin % 

From the first expression, one finds that after integrating by parts that 

C{a) = + ^a* + O (a^ 
as a — > 0. The second expression gives the limiting behavior 

C(a) =a + 0(a- 1 ) 

as a — > 00. 

We make some brief remarks on notation: in this paper we use an integral with a subscript in parentheses 
to denote the contour integral in the complex plane with positive orientation along the line Re(s) = c, 
and use the standard notations < s , o x {-) and O x {-), where the subscript denotes that the implied constants 
depend on the parameter x and are otherwise absolute. We have made an effort to be explicit about the 
dependence of our error terms on the various parameters involved in this paper. 

Finally, we would like to thank the number theory community at Stanford for the stimulating academic 
environment, in particular, Robert Rhoades, Xiannan Li, Bob Hough, David Slier, Akshay Venkatesh and 
above all, Kannan Soundararajan, discussions with whom influenced the final form of this paper. 



f. Connections between the results via Eisenstein and Multiple Dirichlet Series 

There are many parallels between our work and that of Conrey, Farmer and Soundararajan, and moreover 
both results can be interpreted as averages of Fourier coefficients of Eisenstein series. We first mention some 
interesting similarities. 

In both our work and that of Conrey, Farmer and Soundararajan, the result when one parameter grows 
more rapidly than the other can be deduced by evaluating the long sum first and then making trivial estimates 
on the short sum. Indeed, if Y 2 is very small compared to X, Sarnak's solution to the shifted convolution 
problem (|2|) yields 

h<Y n<X 

In the case of the Jacobi symbol, one can use the Polya- Vinogradov inequality to show that 



|f^ + O e (^) ifm = D 

O I 7712 log m ) if 771 7^ □ 



n<X 
n odd 



(-) 



m<Y 
i=odd □ 



from which it follows that 

S(X,Y) = 



,T 2 ' 

and similarly if the roles of n, m are reversed, see pQ. The case of shifted convolution sums when Y 2 is much 
larger than X follows even more simply. Indeed, we have that 



^rXY^ +O e (ytlogF + Y^X £ + XlogY 



n<X 



E E ^n)X f (n + h) = -^^lx + O f (xl 



See for example, [BJ. Then we find 

linU,^ -I- M = -£i 

2C(2) 

by squaring and considering the off-diagonal terms. Observe that these easily-derived estimates correspond 
to the limiting behavior of the transition regions in cither of the more general theorems stated above. 

At the transition phase (i.e. on the blowup) the asymptotic behavior in both theorems depends crucially on 
the automorphic nature of the modular form / or the Jacobi symbol. Conrey, Farmer and Soundararajan's 
proof is based on the Poisson summation formula, see section 3 of their paper [T], which is the essential 
ingredient in proving the functional equations of the Dirichlct L-functions L(s,x), i.e. GL\ automorphy. 
Likewise, one might consider Lemma 1 of section 3 of this paper to be the essential step which brings out 
the behavior of shifted convolution sums on the blowup. The crucial step in the proof of Lemma 1 is the 
application of the functional equation of the L-function L(s, f x /), which follows from the automorphy of 

The Theorem proved by Conrey, Farmer and Soundararajan has one very surprising feature which is 
not apparent in our work. The function C(a) appearing in their result is once continuously diffcrcntiable 
everywhere, but it is twice diffcrcntiable at aQ if and only if a — 2p/q with p and q both odd, see section 
6 of their paper. It is necessary that the sum S(X, Y) has a sharp cut-off for C(a) to have such strange 
differentiability properties. In our work, we do not prove an asymptotic result for a sharp cut-off, but we do 
have some flexibility in our choice of cut off functions. In section 3, we see that the asymptotic we obtain 
on the blowup has a main term which depends on the particular shape of the cut-off function we choose. 
This is unusual in analytic number theory. It remains possible that averages of shifted convolution sums 
in the shift aspect when considered with a sharp cut-off also have strange differentiability properties at the 
transition region. 

The most essential connection, however, is that both theorems arise from averages of Fourier coefficients 
of Eiscnstcin series, and the transition regions can be understood by studying the resulting multiple Dirichlct 
series. Wc first discuss the case of characters. We have that 



1 f f X s V w 

S ( X ' Y ) = 7^2 / Z ( s ' w ) dsdw 

(2in) 2 J, c) J, c) s w 



where c > 1 and 

Z(s, w) 



'(c) J (c) 

(tt) 



E ■ 

m,n>l 
m,n odd 



This is perhaps the first example of a multiple Dirichlct series, and information about the analytic proper- 
ties of Z(s,w) would determine the asymptotic behavior of S(X,Y). Goldfcld and Hoffstcin in [3] derive 
the analytic properties of Z(s,w), which crucially follow from studying weight ^ Eisenstein series for the 
congruence subgroup To (4). More precisely, they study the modified multiple Dirichlet series 

Z+{s,w)= w — ' 

m>l 
m squarefree 

where 

' m = I (mod 4) 



Xm(n) := 



m = 2, 3 (mod 4) 



and find that it has poles along w = 1 and w = | — s. In fact, if Ei(z, s) is the weight ^ Eisenstein series 
at the cusp (see Goldfeld and Hoffstein for details), it has the Fourier expansion 



Ei(z,s) = E a m {s,y)e(nx) 



m>l 



where 



is essentially L(2s,Xm) times a i^-Bessel function. One finds that 

a m (s, l/X) ds 



E - 

^ 2m 

m<Y 
m odd 



(I) 



evaluates to a double mellin inverse of Z(s, w) with some additional factors. Ignoring covergence issues, or 
after sufficient smoothing, one can see how the transition region arises directly from Z(s,w). For simplicity, 
we work with Z+(s,w). We have 



E 



E 

m<Y n<X 
m squarefree 



1 f f , X s Y l 



Xm(n) = , ., 2 / / Z+(s,w) dsdw 



(c)J(c) 



1 f , Y w , , s Xi (Y 

— — / c(w)X he, (10)70 — — 

2mJ (c) y ' w +y y (| - w )w \X 



1 f f X s Y w 

+ / ,, / / dsdu; 

7(c) y(|-Re(to)- e ) s w 

where c(w) and c+(w) are the meromorphic functions arising in Theorem 1 of [4]. They have poles at w — h, 
with res^^i c(w) = — res t0= i c+(w) because of the intersection of the two singular divisors at (s, w) = (|, 1). 
Hence, if Y/X — > 00, the terms involving c(w) and c+(w) cancel, and if Y/X — > 0, the term involving c(w) 
becomes a main term, and the one involving c* + (w) becomes an error term of size 0(Yi). UY/X remains 
constant, the term involving c*j_(w) gives the behavior on the blow-up. 

Our result on the following pages can be seen to arise from an identical situation. Sums of the form (fTJ) 
are essentially an average of Fourier coefficients of |/| 2 . Our method to treat these sums is to take a spectral 
expansion into Eisenstein series and Maass cusp forms, so one is led to compute a sum similar to 



E- 



b m (s, l/X) ds, 

(i) 



m<Y 

where 

E(z,s):= 2J Im(7z) s = ^ &„(s, y)e(nx) 

is the standard weight real analytic Eisenstein series. The integral along Re(s) = | here arises from 
averaging over the continuous spectrum in our decomposition. In this context, the analougc of Z(s, w) is 



EE 



ib—h 



IT 2 ^ = C ( W - S+1 2) (: ( W + S - 1 2 



While not a multiple Dirichlet series, it is still a zeta function in two variables and has singular divisors at 
w = s + I and w = | — s, the intersection of which gives rise to a transition behavior exactly as discussed 
above. This will be carried out in explicit detail in the remainder of the paper. 

Finally, one should note the connection between Ricmann's nondiffcrcntiable function and thcta functions 
on the real axis, see the interesting paper of Duistcrmaat [2]. 



2. Preliminaries 

Our main approach to the Main Theorem is to take the Petersson inner product of y k \f\ 2 against an 
incomplete Poincare series P/ l (-| , 0). This approach was first introduced by Selberg [15j . and has been suc- 
cessfully used by many other authors to study shifted convolution sums in the past (for an overview, see 
|13j). Throughout this paper we set T = SL^^Z), and let H denote the upper half plane with its hyperbolic 
metric. We work in the Hilbert space L 2 (T\H) of square integrable measurable functions with the Petersson 
inner product 

(u, v) = / u(z)v(z) df+z. 
Jr\-H 
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The symmetric operator 

A = - 2 (— + — 

^ \ dx 2 dy 2 

acts on the subspace of smooth functions and moreover has a unique self-adjoint extension to all L 2 (T\H), 
see Iwaniec [11] chapter 4. Given a classical holomorphic normalized cuspidal eigenform / of weight k, set 
F(z) := f(z). We have that |F| G L 2 {T\H), but on the other hand it is no longer holomorphic. Let 4>(y) 
be an infinitely differentiable compactly supported function on M>o, and let denote the stabilizer in T 
of the cusp at infinity. Then we define the incomplete Poincarc scries 

P h (z\ip) := Y e(hjz)ip(lm(jz)) , 
7er oc \r 

which is a smooth and bounded function on T\H. By unfolding the inner product 



(FP h (-\i>),F)= / y k \f(z)\ 2 P h {z\i>)d^z 

JV\H 

on the Poincare series, we find that 

OO n OO 

(FP h (-\ij),F) = y2x f (n)X f (n + h)(n(n + h)) ! ^ 1 / ^y)e^ n+h V" 2 dy, 
i Jo 

n—l u 

i.e. this inner product is a smoothed shifted convolution sum with cut-off function given in terms of an 
integral transform (similar to the Laplace transform) of ip. The behavior of ip{_y) as y tends to is crucial 
to control the length of the shifted convolution sum. In connection with the previous section, it should be 
noted that if if) were a delta function, then taking the inner product against Ph(z\ip) is equivalent to taking 
the h-th Fourier coefficient. 

Let Uj be a complete orthonormal system of cusp forms which are eigenfunctions of the Laplace operator 
and all Hecke operators. Because we are only working in level 1, we need not worry about old forms or the 
Hccke operators whose index divides the level. Define the real analytic Eiscnstcin series by 

E(z,s):= Y, Im ^ z ) S 
7er oc \r 

for Re(s) > 1, and in general by analytic continuation. For each s ^ 0,1, the Eisenstein series are also 
eigenfunctions for the Laplace operator and all the Hecke operators. We have then that 

Y(P h (-W),u j )u j (z) + — / (P h {-\iP),E(;l/2 + it))E(z,l/2 + it)dt 

j=l J -oo 

converges to Ph{z\ip) in the norm topology on L 2 (T\H), see Theorems 4.7 and 7.3 in [TT]. Then we have 
that 

(3) (FP h (-W),F)=Y( p h(-\i>),u j )(Fu j> F) + — / (P h (-\i>),E(;l/2 + it))(FE(;l/2 + it),F)dt, 
i=i J—oo 

and we will see later that the convergence is absolute and uniform in h. 

If u G L 2 (T\H) and Au = Xu with A = s(l — s), then u(z) has a Fourier expansion given in terms of the 
A'-Bcssel function K v (z), which is an exponentially decaying solution to the differential equation 

z*r + zf>-(z 2 + v 2 )f = 0. 

We will primarily be interested in the X-Besscl function for purely imaginary v, and note some useful 
properties of these functions: first, that K v (z) is real for real z, second, that as Im(z^) — > oo in a fixed 
vertical strip, K v (z) is decaying exponentially, and last, that for v = it, with Kn{z) has a branch cut, 

which we take to be along the negative real axis in the z-plane. As z — > 0, we have that 

sinolog z/2) 

\Kit{z)\ 

so long as one avoids the branch cut. 



r(l + it) sinh(Trf) 



3. Proof of Theorem 

We now proceed to the proof of the Main Theorem of the paper by summing the right side of © over h. 
Pointwise, the largest term comes from the discrete spectrum (see [2]), however, on average, the continuous 
spectrum dominates. We start with the continuous spectrum contribution to ([3]). 

3.1. Eisenstein Series. The Eisenstein series E(z,s) has a Fourier expansion given by 

<p(n,s)W s (nz), 



where 



W s (z) = 2yiK s _±(2TT V )e{x) 



If n = and s/:| then 
and if n ^ 0, the n-th Fourier coefficient of E(z, s) is given by 



V (n, S ) = «2»)-'|„r 1/2 £ (f 



s-1/2 



where £(s) = n s / 2 T(s/2)((s) denotes the completed Ricmann zeta function which has the functional equa- 
tion £(s) = £(1 — s). If s = h + it, we find by unfolding the Poincare series that the inner product 

WW,^,l/Hit)) = ^E(J) S4 f^±e-^yK s _,(2«hy)dy 

where we have used that £(2s) = £(2s) = £(2 — 2s) = £(2s — 1). We can also unfold the second inner product 
on the Eisenstein series. Following Iwaniec [10] chapter 13, set 

L(s, f x /) := C(s)L(s, sym 2 /) = C(2s)L( S , f ® /) = C(2,s) V ^T^> 

n>l 

where the last equality is valid only for Re(s) > 1. This L-function admits the functional equation 

A( s ,/x/)=A(l-s,/x/), 

where 

A(s,fxf) = L 00 (s,fxf)L(sJxf) 

with 

Loo(s, / x /) := (2^)- 2s r( s )r( s + k - 1). 

By unfolding when Rc(s) > 1, and in general after analytic continuation, 

Going back to the spectral expansion (|3|) and pulling these two inner products together, we have that the 
Eisenstein series contribution to (F Ph(-\ip) , F) is 

fMY) (4^-i 2 ™y (|) £(2 S )£(2 S -l) a ^Uy J j/3/2 s 2 K y> 

The completed L and zeta functions have the same number of gamma factors in the numerator and denom- 
inator, and Bcssel function decays rapidly as |Im(s)| — > oo, so the contour integral converges rapidly. It is 
interesting to note that the s integral only makes sense because £(s) -1 has no poles and at most polynomial 
growth on the Re(s) = 1 line, i.e. due to the prime number theorem. Indeed, one sees that the contour 
is constrained between the poles in the critical strips of the two £(s) -1 functions appearing here. Due to 
the mysterious nature of the residues at these poles, shifting contours seems to be a futile approach to 
understanding the asymptotic size of Ef^ty). 

7 



Introducing a sum over h clears this obstruction. We have to compute 

h<y v ; *' u » "^2) sv ;sv ' h<Y ab=h 

and can evaluate the sum over h by standard techniques. We have by adapting Theorem 12.4 from |17j (due 
to Van Der Corput) that 



^ll + l Re(s)— §|+e 



EE© -cm^+cp-^ + o. 

h<a; ab—h z ^ 

where the implied constants depend at most polynomially on \s\. It should be noted that the error term 
here is not the best currently known, however, it is sufficiently small as to not contribute to the final result 
of this paper. By partial summation, 



h<Y ab=h 
i-Y 



= C(2s) / u s -h- 27Tyu K s _i(2Tryu) du + ({2 - 2s) [ u^~ s e~^ yu K s _i {2-nyu) du 
Jo 2 Jo 2 

+ O s . e (e- 2nyY K s _i_{2iiyY)Y& + \ Kc ^-^ + ^ . 

Note that the two integrals appearing in the displayed equation are interchanged under the transformation 
s < — > 1 — s, by symmetry of the the Besscl function. The contour integral over s in ^2 Ef,h{$) is also 
symmetric under s < — > 1 — s so that these two integrals are identical in the overall sum, and we need only 
work one of them out. The first integral can be evaluated explicitly, and the answer will be in terms of 
hypergeometric functions. The confluent hypergeometric function that will appear below is defined by the 
power series 



iFi(a,b,z) = ^2 



a (n) z n 



„(")n! 

n— 

where 

a 



(") = a(a + l)(a + 2) • • • (a + n - 1) = r( ° + . n) 

r(a) 

is called either the 'rising factorial' or 'Pochhammer symbol'. The theory of hypergeometric functions is 
developed in detail in [3]. It is entire on C separately in each variable except for simple poles at b = 
0, — 1, —2, ... by the absolute and uniform convergence of the defining series, hence it is meromorphic on C 3 . 
We have that the residues at these poles are given by 

res l F l (a, b, z) = + t^" 1 ]" , (a + n + 1, n + 2, z), 

b=—n I (a)I {n + 2)1 [n + 1) 

see Gradshteyn and Ryzhik [5], 9.214. 
We have the formulae 

T . , \ 7r i" 'J- V {iz) - i~ v J u (iz) 
I sin -kv 

e- z T ( S+ j) J s _i(iz) = i..Fi(s,2s,-2z), 

/ iz_ \ s 2 2 

2 ; 

^-1^1(0, b, z) = jiF^a + 1,6+1, z), 
dz 

and 

j- {z b -\F 1 {a, b, z)) =(b- l)z b ~\F 1 (a, b - 1, z), 

where J„(x) is the J-Bessel function, the second formula can be found in [S], 9.215 #3 and the last two 
formulae can be found in [5] section 2.1.2. From these it follows that 

u s-± e -2. V u K {2lTyu) du = ( Af {y, s) + B f (y, Y, s) + C f (y, Y,s)), 

A{iry)2s 



where 

A f = A f {y,s) = (iry)- s T(l/2 + s) 
B f = B f (y, Y, s) = (^)- s r(l/2 + a^F^-a, 1 - 2s, -4nyY) 

and 

C f = C f (y, Y, a) = (TryY 2 ) s T(l/2 - s^F^s, 1 + 2s, -AnyY). 

The C(2s) from the evaluation of the sum over h cancels against the £(2s) in the denominator of (QJ, 
eliminating its poles. Given that Af + B j + Cf is holomorphic in s past the Re(s) = line, we are free to 
pass the contour to the left. Picking up a residue at s = we get 

v-p m r(fc) X(l,sym 2 /) f°°i>(y) A 

^ EfM = 2C(2) L ~^ dV 



h<Y 



1 f 1 f A(sjx f)(A f (y,s)+B f (y,Y,s) + C f (y 2 Y 1 s)) ds 



(4ir) k -i Jo V 2 ^iJ(-a) f(2s-l) 7r-T(a + l) 



V Jo y 2 

where < a < ^ and we have again made use of the prime number theorem in estimating the error term. 

The integrals over y converge at oo, so the asymptotic size of FjmW depends only on the behavior of 
ip(y) for small y. We now make some estimates assuming that y is small. The remaining contour integral 

1_ f A(s,fx f) (A f (y,s) + B f (y,Y,s) + C f (y,Y,s)) ^ 



2tu ./(_„) £(2a-l) 7r-r(*+l) 

is a sum of three terms coming from At, Bf and Cf. The term coming from Af is -Cfc y^ as y —> 0. 
The hypcrgcomctric function appearing in Bf is bounded above and below by universal constants in the 
half-plane Rc(s) < and when < yY < 4ir, thus term coming from Bf is also <C J/ 2 , uniformly in Y. Thus 
it remains to inspect the term coming from C / . Explicitly, let 

2W(-a) £(2s-l) 7T T(s + 1) 



Thus 



j? m r(fc) L(l,sym 2 /) f 00 ^(y) 1 , vW 



In the transition region, Cf_i(y, Y) is the crucial term. 

Lemma 1. Suppose that y,Y £ R>o with Y becoming large and y becoming small. Then 



Cf,i(y,Y) = ^j=cf (4^y 2 ) + o k (yi) 



where 



and 



c/(a) == \ol Xf(n) 2 W k {n 2 na), 

n>l 



i f r( s + fe-i)r( s -i) _ s 



2™ 7ri+a) r (2 - s) 



for any fixed a > . 



Proof. We apply the functional equation for the L-function to find that 

9 9 1 f L(s,f X /) T(s + k- l)T(s- i) , w , 

C Ll {y,Y)=* 2 yY 2 — , ~ ^1^(1 - S , 3 - 2 S , -4rr y r) 4^yF 2 da. 

2m J (1+a) C{2s) T(2 - s) 

From the definition one finds that i-Fi(l — s,3 — 2s, u) = 1 + O s (u), which we use to eliminate the hy- 
pergeometric function from the above expression. We proceed in two slightly different ways depending on 
whether yY 2 becomes large or becomes small. If yY 2 remains bounded, either approach is acceptable. First, 
assume that yY 2 is becoming small. In this case, choose a = j, and observe that the s-dependence in the 
hypergeometric function is uniformly bounded along the line Re(s) = |. Together with the rapid decay of 



the integrand of Cf,i(y,Y), this gives us that 

<fc vHyY 2 )' 



CfAy ' Y) ~2~^ Cf 



If yY 2 becomes large shift the line of integration to the right, past the pole of the hypergeometric function 
at s = | to a = |. The contribution to C/.i(y, Y) coming from this residue is <Cfc y 5 1 uniformly in Y, and 
the s-dependence in the hypergeometric function is uniformly bounded along the line Re(s) = |. We find in 
this case that 

C fAv,Y)-^=c f (AnyY 2 ) 



In either case, we obtain the error term stated in the Lemma. 



□ 



Lemma 2. The function c/(a) defined above is C°°(R>o)- As a — > oo it decays faster than any polynomial, 
and as a — > 

r(fc)L(l,sym 2 /) 

c/(a) = — 2cm — + Ek[a) 

where Ek(a) = o k unconditionally, and Ek(a) = Ok, e |a^~ E j for any e > assuming the Riemann 

hypothesis for the classical Riemann zeta function. 

Proof. The Wk(x) defined above is C°°(R) and its integrand has no poles to the right, thus Wk(x) is rapidly 
decaying as x — > +oo. After differentiating the series for Cf(a) in the second line of Lemma 1 arbitrarily 
many times, the resulting series for (a) converges absolutely for any a > 0, so Cf <G C°°(R>o). The rapid 
decay of c/(a) follows from that of Wk(x). By shifting the line of integration in the definition of c/(a) to 
the left, we investigate the behavior of c/(a) near a = 0. The main term comes from the residue of the pole 
of L(s, f x /) at s = 1, and the error term is estimated by pushing the contour just past the Re(s) = | line, 
or to the Re(s) = \ + e line if one assumes the Riemann hypothesis. □ 

Thus we have 

Proposition 1. For Efh(i>) and c/(a) as defined above, we have that 

The term involving || is smaller than the remainder terms coming from Maass forms, as we will see in 
the next section. 

3.2. Maass Forms. The discrete spectrum of A is spanned by Maass cusp forms. The Hccke algebra acting 
on L 2 (T\H) is defined to be the algebra generated by the commuting self-adjoint bounded operators T„, 
where for u S i 2 (r\H), define T n by 



(MM:- £E £ .(=±i) 

V nr l— „ l. I J <7\ \ / 



ad=n b (mod d) 
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These operators commute with A as well, so in fact our basis of Maass forms can be taken to be eigenfunctions 
of the Hecke algebra as well, and we denote the Hecke eigenvalues of the Maass form Uj by X Uj (n) . A Maass 
form of Laplace eigenvalue Xj = Sj(l — Sj) = j + is cuspidal, so it has a Fourier expansion of the form 

u j( z ) = E a Uj (n)W Sj (nz), 

where 

W Sj (z)=2yiK itj (2Try)e(x). 
For T = SL2(Z), it was known to Sclberg in the early 50s that the smallest Laplace eigenvalue Ai is 
> j, hence tj € K. For a proof of this fact, see [7], chapter 11. Computationally, it has been verified that 
t\ = 9.53369526 . . ., see for example [8]. To apply the spectral theorem wc must assume the normalization 
\\uj 1 1 £ 2 = 1, in which case the Fourier coefficient and Hecke eigenvalue are related by 

i_ 

/ COSllTTt, \ 2 

\ 2|n|L(l, synruj) / 
where the symmetric square L-function appearing here is defined 

The Hecke eigenvalues conform to the bound \X Uj (n)| < d(n)n , where the generalized Ramanujan conjecture 
implies that 9 = is admissible, and the best known unconditional bound is due to Kim and Sarnak [12j . 
which gives 6 = 

We call the Maass form contribution to ([3]) 

oo 

M M (V0 :=^(P,(#), % -)(F % -,f). 
i=i 

By unfolding we have 

<^(#)>^> = ( rn^ O '^W r^e-^K itj {2^hy)dy 
\L{l,sym 2 Uj)J J V 2 



so that 



E M fM = r^£( rn 0Sh 1 3 y ^ F ) E K 3 {h)^ hy K^hy)dy 



7/2 ^-^ V L(l, svm 2 w 7 ) , 

h<Y " u » 3 = 1 v v ' J 3// h<Y 

Lemma 3. The spectral sum 

appearing in ^ Mf_h(?p) converges absolutely. 

Proof. There are three factors in the summand: that which involves cosh7rfj , the inner product, and the 
sum over h. Wc show that the first of these two balance each other, and then show that the sum over h 
decays rapidly in \tj\, uniformly in the other variables. To study \(Fuj,F)\ we will use a beautiful formula 
of Watson [18], but follow a classical work-out of it from Soundararajan's paper [16] ■ In that paper both / 
and Uj are normalized to have mass 1, but in this paper we take / to be Hecke normalized so that Watson's 
formula is 

UFv FM2 _ i / r(fc) y i(r\?m 2 A(i,/x/x % ) 

K J ' 8\(47r)K C(2) J L^Lsym^ACLsym^)' 

where 

A(s, / x / x Uj) = ioo(s, / x / x Uj)L{s, f x / x Uj), 

li 



L^s, sym 2 /) = tt 2 T 



s + 1 



s + k - 1 



s + k 



A(s,sym 2 Uj) = L 00 (s,sym 2 Uj)L(s,sy"m 2 Uj), 



ioo(s,/x / XMj ) = ^- 4s []r 



s + — 1 ± ii. 



s + /c ± it, 



s + 1 ± it. 



s ± it. 



and 



L(s,fx f x Uj) 



Loc^symS) - -- 3s/2 r (^) r © r (i±a 

tjA «/(rt 2 ^(p) \ _1 A _ a Uj (p) \ ~ 2 L _ l3 } {p) 2 a Uo {p) \ _1 



V P s J V P s J V P s / 



The archimcdian parts evaluate to 

ioo(f,/ x / x 



47T' 



|r(fc - ± + itj 



L 00 (l, S ym 2 /) 2 ioo(l,sym 2 ^) r(fc) 2 
after repeated application of the duplication formula. Using this, Watson's formula simplifies to 



(5) 



\(F Uj ,F)\ = V2 l - 



|r(fc- |+^)| fLfofxfxu,) 



(4tt) a 



L(l, sym 2 itj 



for / Hccke normalized, and Uj mass 1 normalized. By applying Stirling's formula, we find that \T(k — \ + 
iij)|(cosh7rtj)5 is polynomially bounded as \tj \ — > oo. Together with standard convexity bounds in the \tj\ 
aspect for L(^,f x / x uj), we find that 



cosh irt,; 



L(l, sym 2 z 



\(Fu 3 ,F)\ 



is polynomially bounded as \tj\ gets large. 

Now we turn to the sum over h. We have by a "folklore" result written down by Hafner and Ivic [5] that 

h<Y 

where 9 = or = ^ as above, and where the implied constants depend at most polynomially on \tj\. As 
on the Eisenstein series side, the conjectural truth is O u ^(Y^ +E ), but this seems very difficult. By partial 
summation 



Y K{h)e~ 2 ^yK ttj {2^hy) <„. e~ 2 ^ Y K itj {2i:yY)Y + £ u* 

h<Y "2 



(1+9) 



^-e~ 2 ^K ltj (2nyu) 



du. 



where the implied constants again depend at most polynomially on \tj\. We have that \K it (u) | ~ it | p ^ " ^ )°f ^ ( x t ) 
for small u, thus after taking derivatives, changing variables, and using the power series expansion for the 
lower incomplete gamma function, we have that 

(6) Y, Uh)e- 2 ^K ltj (27rhy) < P(t 3 )e-^Y^ 1+6 \ 

h<Y 

uniformly in y, where P(t) is a real- valued function on K that grows at most polynomially as |t| becomes 
large. From these estimates together with Weyl's law (see, e.g. Iwaniec QT]) 

Y 1 = V0l ^ ) T 2 + Q £ ((l+T) 1 + £ ) 



the Lemma follows. 



□ 



We apply trivial estimates along with ((6]) to obtain 
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Proposition 2. For Mf^fjp) defined above we have 



h<Y J ° 



Mv)\ 

3 

V 2 



dy. 



The application of trivial bounds is justified by the absolute convergence given by Lemma 3. 



4. The Main Theorem and its Corollaries 

Drawing together the propositions from the two preceding sections, we obtain 

Main Theorem. Let ip be any measurable function on R>o such that the incomplete Poincare series 

Ph(z\tp) := ^2 e(hjz)^(lm('yz)) 
7er 00 \r 

is a smooth and bounded L 2 function on T\H. Denote the shifted convolution sum 



S f (i>,Y) = J2J2x f (n)X f (n + h)(n(n + h)Y 
Then we have that 



h<Y n>\ 



^y) e -Mn+h)y y k-2 dy _ 



T(k)L(l,sym 2 f )\ ^(y) 
2C(2) 



■dy + O k [yia+«) 



Mv)\ . 

— dy 



where Cf(a) is the function defined in Lemma 1 and whose properties are given in Lemma 2, and 6 = or 
depending on whether we assume the truth of the generalized Ramanujan conjecture or not. 

It should be noted that the size of the exponent of Y in the error term here depends on the sharp cut-off 
in h, and could be made smaller if we were to smooth that sum. Now we make some choices for ijj, and 
record the results as corollaries. First, we give the result stated in the introduction. 

Proof of Corollary 1. Let tp be a smooth approximation to a point mass. Specifically, let ip be smooth, 
non- negative, supported on a set of radius X~ A about the point y = and have mass 1. Then for any 
continuously differentiable function <f> on K>o, we have 



/ 1 



i>(y)<l>(v) dy 







< 


4 1 ) 





x- 



where the implied constants are absolute. First, let <p(u) = u k 2 e 47r ("+' I ) M j so that \<j>' (4^) | <^fc (n + 
h)X 3 ~ k e~~^~ , and thus we find that 



(n(n + h)) 2 _n±h 
e x 



h<Yn>l 



X 



k-2 



hence the difference between the left hand sides of the Main Theorem and Corollary 1 is <C £ X e . Secondly, 
let 

r(£0L(l,sym 2 /A 1 



0= (c^uY 2 )- 



2C(2) 



From the definition of c/(a) as a sum one sees that |cy(a)| <C/ oT l as a — > 0, so that yjjf "C/ X 3 . 

The error term can be treated similarly. Hence, the difference between the right hand side of the Main 
Theorem and the right hand side of Corollary 1 is -C/ X~ 2 with this choice of tp. 



□ 



Let 



denote the incomplete gamma function. 



T(s,x) 
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Corollary 2. Let 



Then 



i\ ) i\ \ n + hl r(fc-i) 

h<Yn>l v 7 v ; 



L(l,sym 2 /) F 2 p ^(u) f^|(i+e) 

where the integral appearing here either cancels the main term if — approaches 0, or decays rapidly as a 
function of if this parameter grows without bound. 

Proof. Let 1 >X denote the indicator function of the open set (x, oo) C R. In similar fashion to the previous 
proof, let ij)(y) be a smooth approximation to 

and compute the answer. □ 

The spectral theorem as used in this paper holds for L 2 functions which are also C°° and bounded, 
however, it should also hold for a much wider class of functions, and thus our Main Theorem should in fact 
hold for a much wider class of functions. In practice however, one may always obtain a corollary for a specific 
choice of ip by elementary arguments similar to the proof of Corollary 1, so we do not pursue the problem 
of expanding the class of functions for which the Main Theorem holds. 
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